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During recent years our knowledge of the thermodynamics of strongly interacting el- 
ementary particles greatly advanced. Lattice calculations now allow to study also the 
thermodynamics at non-zero quark chemical potential {jl q ). The different approaches 
developed for this purpose 111 0. S 0. S are still limited to the regime of high temper- 
ature and small values of the chemical potential, T > 0.9T C , }lq/T < 1. They, however, 
allow already to analyze the density dependence of the equation of state in a regime rel- 
evant for a wide range of energies explored by heavy ion experiments and may even be 
sufficient to establish or rule out the existence of a second order phase transition point 
in the QCD phase diagram. The existence of such a critical point as endpoint of a line 
of first order phase transition that separates at low temperature the low and high density 
regions in the QCD phase diagram, is suggested by many phenomenological models. 
For small values of H- q /T lattice calculations suggest that the transition from low to high 
temperature is not a phase transition; the transition during which bulk thermodynamic 
quantities, e.g. the energy density, change rapidly and the chiral condensate drops sud- 
denly, is a continuous, rapid crossover. It thus has been speculated [6] that a 2 nd order 
phase transition point exists somewhere in the interior of the QCD phase diagram. 

The generic form of the QCD phase diagram is shown in Fig.[U Although this phase 
diagram is well motivated by model calculations little is known quantitatively about 
it from lattice calculations. All calculations performed so far for QCD with two light 
quarks with or without the inclusion of a heavier strange quark find a smooth but rapid 
crossover from the low to high temperature regime. However, these calculations also 
did not provide evidence for the expected universal critical behavior in the light quark 
mass limit, which for 2-flavor QCD should be that of a 3-dimensional 0(A) or 0(2) 
spin model. It generally is expected that calculations closer to the continuum limit are 
needed to unreveal these universal features of the QCD transition. In view of this missing 
evidence for O(N) scaling, it has been argued recently that the transition in 2-flavor QCD 
could also be a weak first order transition [7], in which case the entire transition line in 
the \l-T phase diagram could be a line of first order transitions. 

We will start our survey of lattice calculations on QCD thermodynamics in the next 
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FIGURE 1. Sketch of the QCD phase diagram. 



Section by discussing recent studies of the QCD equation of state at vanishing and non- 
vanishing chemical potential. Section 3 is devoted to a discussion of the status of studies 
of the 2 nd order phase transition point (chiral critical point) in the QCD phase diagram. 



Most information on the structure of the high temperature phase of QCD and the 
nature of the transition itself has been obtained through lattice calculations performed 
in the limit of vanishing baryon number density or vanishing quark chemical potential 
(}l q = 0). This limit is most relevant for our understanding of the evolution of the early 
universe and also is the regime which can be studied experimentally in heavy ion colli- 
sions at RHIC (BNL) and soon also at the LHC (CERN). The experimental accessibility 
of this regime of dense matter also drives the desire to go beyond a qualitative analysis 
of the QCD phase transition and to aim at a numerically accurate determination of basic 
parameters that characterize the thermodynamics of dense matter at high temperature. 

The transition to the high temperature phase of QCD is related to the restoration of 
chiral symmetry as well as deconfinement; the vanishing of the chiral condensate at 
the transition temperature Tq and the sudden liberation of quark and gluon degrees of 
freedom is clearly visible in Fig. As can be seen the temperature and quark mass 
dependence of the light quark chiral condensate in QCD with two light and a heavier 
strange quark is similar to that of QCD with 3 degenerate (light) quarks. In the high 
temperature phase, on the other hand, bulk thermodynamic observables, e.g. the pressure 
shown in Fig.|2left), clearly reflect the number of light degrees of freedom and also are 
sensitive to the heavier strange quark mass 1 . 



1 In the calculation of the pressure of (2+l)-flavor QCD shown in Fig. fright) the strange quark mass 
has actually been taken to be proportional to the temperature, i.e. m s /T =const. If one keeps instead m s 
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FIGURE 2. The light quark chiral condensate in QCD with 2 light up, down and a heavier strange 
quark mass (open symbols) and in 3-flavor QCD with degenerate quark masses (full symbols) The 
right hand part of the figure shows the pressure calculated in QCD with different number of flavors as well 
as in a pure gauge theory [9]. Note that (2+l)-fiavor QCD here refers to QCD with two light quarks and a 
heavier (strange) quark with a mass proportional to the temperature, m s ~ T. 



Recent studies of the equation of state concentrated on calculations with an almost 
realistic quark mass spectrum performed along lines of constant physics, i.e. with light 
and strange quark masses fixed in units of hadron masses rather than the temperature 
Till . To reduce cut-off effects induced by the finite lattice spacing larger temporal 
lattices (N t = 6) and so-called fat links have been used. At present these calculations 
still have been performed with rather small spatial volumes, V l ^T ~ 2. Close to Tq, 
where the correlation length becomes large, a more detailed analysis of the approach to 
the thermodynamic limit thus still has to be performed in the future. 

Partly these calculations still have been performed with unimproved fermion actions 
iflflll which makes it difficult to control cut-off effects and draw quantitative conclusions 
on the approach to the ideal gas limit at high temperature. Nonetheless, these calcula- 
tions support earlier findings on the temperature dependence of the pressure and energy 
density in the transition region. In particular, they show that the presence of strange 
quarks has little influence on the thermodynamics in the vicinity of Tq. In Fig.|3]we com- 
pare the recent calculation of the pressure in (2+l)-flavor QCD II 1 Qfl and earlier results 
for 2-flavor QCD II12L 1 1 3fl which both have been performed with unimproved staggered 
fermions. The good agreement in the vicinity of Tq suggests that the strange quark con- 
tribution to the pressure is small. Only for T > 1.57b differences show up; the positive 
strange quark contribution to the pressure in (2 + 1) -flavor QCD becomes sizeable 2 . 

The good agreement between recent calculations in (2+l)-flavor QCD, performed 
with smaller quark masses and smaller lattice spacings iTToL [Till , and earlier results 
§01 EI is quite reassuring. These calculations confirm that thermodynamics in the 
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fixed to its physical value the strange quark effectively becomes light in units of the temperature and the 
pressure will gradually approach the high temperature limit of 3-flavor QCD. 

2 As the relative size of the strange quark contribution to thermodynamic quantities changes with temper- 
ature it is obvious that cut-off effects can not be taken into account by rescaling the pressure or energy 
density with a universal, temperature independent factor fPTOll . This underlines the need of improved ac- 
tions for thermodynamic calculations. 




FIGURE 3. Cut-off dependence of the pressure calculated with the standard staggered fermion action 
on lattices with temporal extent N r = 4 and 6 in 2-flavor QCD (T3IH] and (2+l)-flavor QCD (kJ (left). 
The right hand figure shows energy density, pressure and / — e — 3p for (2+l)-fiavor QCD calculated on 
lattices with temporal extent jV t = 4 and 6 with an &(a 2 ) improved staggered fermion action full . 



high temperature phase is rather insensitive to changes of the quark mass; a reduction of 
the light quark masses by almost an order of magnitude does not lead to drastic changes 
in the energy density at the transition point and in the high temperature phase. Moreover, 
they show that the transition itself is not strongly influenced by discretization errors, 
which in the staggered fermion formulation show up prominently in the distortion of 
the Goldstone modes; reducing m q and thus the masses of light hadrons as well as 
reducing flavor symmetry breaking effects drastically [ 10] does not significantly change 
the energy density at the transition temperature. The estimate, £ c /Tq = 6±2 J2D, is 
consistent with the recent calculations in (2+l)-flavor QCD performed with lighter quark 
masses. The weak dependence of transition parameters on the quark mass reflects the 
importance of numerous heavy resonances that are necessary to build up the particle and 
energy density needed for the transition to occur ifbUl . 

Also the recent studies of the equation of state performed at vanishing quark chemical 
potential suggest that for physical values of the quark masses the transition to the high 
temperature phase of QCD only is a rapid crossover rather than a phase transition 
which on finite lattices would be signaled by metastabilities in bulk thermodynamic 
observables or the chiral condensate. None of the calculations performed so far for QCD 
with two light quarks with or without the inclusion of a heavier strange quark gave direct 
evidence for a first order phase transition. 

The missing guidance from any universal behavior in the vicinity of the transition 
also influences the determination of the transition temperature itself. As the transition 
temperature is determined at quark mass values which are usually larger than those re- 
alized in nature one has to extrapolate to the physical regime. Depending on the ansatz 
used for the quark mass dependence of T c the extrapolations can differ by about 5%. A 
similar systematic uncertainty arises from the calculation of a zero temperature observ- 
able that is used to set the physical units for T c . Thus also the recent determinations of 
the transition temperature illlLll5i1 . which lead to values T c ~ 170 MeV still suffer from 
systematic errors of about 10%. 




Non-zero baryon number density 

Studies of the QCD equation of state have recently been extended to the case of 
non-zero quark chemical potential (/i^). Calculations of bulk thermodynamic quantities 
for fj, q > based on the reweighting approach [16], using the Taylor expansion of 
the partition function H17L 1 1 3L Il9ll . as well as analytic continuation of calculations 
performed with imaginary values of the chemical potential lEoll show that the /in- 
dependent contributions to energy density and pressure are dominated by the leading 
order (fi q /T) 2 correction for parameters relevant for the description of dense matter 
formed at RHIC {}JL q /T ~ 0. 1). Higher order contributions are still only a few percent at 
SPS energies (jj^/T < 0.6). 

We will focus in the following on a discussion of Taylor expansions of the partition 
function of 2-flavor QCD around \i q = 0. At fixed temperature and small values of the 
chemical potential the pressure may be expanded in a Taylor series around \x q = 0, 

w=0 

where the expansion coefficients are given in terms of derivatives of lnZ(r,/i ? ), i.e. 
1 d n \nZ 

c n {T,m q ) = —.-=rj — i — r^- The series is even in (/i q /T) which reflects the invariance of 

Z{T, }i q ) under exchange of particles and anti -particles. The Taylor series for the energy 
density can then be obtained using the relation, (e - 3p)/T 4 = Td(p/T 4 )/dT, 

ji= £ {3c n (T,m q )+c' n (T,m q )) {^ff (2) 

n=0 

with c' n (T, m q ) = Tdc n (T, m q )/dT. A similar relation holds for the entropy density ifl^l . 
The coefficients c n (T, m q ) calculated for a fixed value of the bare quark mass up to n = 6 
are shown in Fig.HJ 

Knowing the dependence of the energy density and the pressure on the quark chemical 
potential one can eliminate \x q in favor of a variable that characterizing the thermody- 
namic boundary conditions for the system under consideration [19]. In the case of dense 
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FIGURE 5. Equation of state of 2-fiavor QCD on lines of constant entropy per baryon number. The 
left hand figure shows three lines of constant S/N b in the QCD phase diagram relevant for the freeze-out 
parameters determined in various heavy ion experiments. The right hand figure shows the equation of 
state on these trajectories using T c = 175 MeV to set the scale. The solid curve in the right hand figure is 
the parametrization of the high temperature part of the equation of state given in Eq.[3] 



matter created in heavy ion collisions this is a combination of entropy and baryon num- 
ber. Both quantities stay constant during the expansion of the system. In Fig.|3]we show 
the resulting isentropic equation of state as function of temperature as well as energy 
density obtained from a 6 th order Taylor expansion of pressure and energy density yj|]. 
The three different entropy and baryon number ratios, S/Nb = 30, 45 and 100, corre- 
spond roughly to isentropic expansions of matter formed at the AGS, SPS and RHIC, 
respectively. It is quite remarkable that p(e) is to a good approximation independent of 
S/Nb', for temperatures T > Tq, or equivalently £ > 0.8 GeV/fm 3 , it is well described by 

P = l(l 

£ 3 V l+0.5£fm 3 /GeV 

which for large £ agrees with a bag equation of state with 5 1 / 4 ~ 260 MeV. 

The insensitivity of the isentropic equation of state on S/Nb also implies that the 
velocity of sound, v| = dp/ d£ is similar along different isentropic expansion trajectories. 
In fact, the parametrization given in Eq.|3]suggests that the velocity of sound approaches 
rather rapidly the ideal gas value, v| = 1/3. In Fig. |6] we summarize results for 
obtained in lattice calculations for a pure SU(3) gauge theory 12 ill , for 2-flavor QCD 
with Wilson fermions 12211 and (2 + l)-flavor QCD with staggered fermions llldl . 




THE CHIRAL CRITICAL POINT 

Various model calculations JfJ suggest that a second order phase transition point (chiral 
critical point) exists in the QCD phase diagram which separates a region of first order 
phase transitions at high baryon number density and low temperatures from a cross-over 
region at low baryon number density and high temperature. Evidence for the existence 
of such a critical point may come from lattice calculations at non-zero quark chemical 
potential by either determining the location of Lee- Yang zeroes IH or by determining 
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FIGURE 6. The velocity of sound in QCD vs. temperature expressed in units of the transition temper- 
ature Tq. Shown are results from calculations with Wilson 1 22] and staggered fermions 1 10] as well as for 
a pure SU(3) gauge theory [21]. Also shown is the resulting deduced from Eq.l3l fl^l . 



the convergence radius of the Taylor series for the logarithm of the partition function 
which directly yields the pressure, p/T = V~ l lnZ aim . 

If there exists a 2 nd order phase transition point in the QCD phase diagram, this could 
be determined from an analysis of the volume dependence of Lee- Yang zeroes of the 
QCD partition function. In any finite volume zeroes of Z(V, T,ll q ) only exist in the 
complex \x q plane with lmfi q ^ 0. Only for V — > °° some of these zeroes may converge 
to the real axis and will then give rise to singularities in thermodynamic quantities. The 
relation between phase transitions and zeroes of the partition function has been exploited 
using a reweighting technique to extend lattice calculations at /l q = to /l q > (H. 
Recent results based on this approach [23] suggest that a critical point indeed exists and 
occurs at Ur = 3jl q ~ 360 MeV. This estimate is about a factor two smaller than earlier 
estimates [ 1] which have been obtained on smaller lattices and with large quark masses. 
This suggests that a detailed analysis of the quark mass and volume dependence lE4ll 
still is needed to gain confidence in the analysis of Lee- Yang zeroes. 

The radius of convergence of the Taylor series is controlled by a singular point in 
the complex \x q plane closest to the origin. It is related to the location of the critical 
point only if this singularity lies on the real axis. A sufficient condition for this is that 
all expansion coefficients in the Taylor series are positive. For temperatures below the 
transition temperature at \x q = this is indeed the case for all expansion coefficients 
calculated so far. The first coefficient, cq, just gives the pressure at \x q = shown in 
Fig.EJnght) and thus is positive for all temperatures. This also is the case for ci, which 
is proportional to the quark number susceptibility at \l q = II25II . 

= IZ'JtY = t dn (t)" ^ = (" + 2 )("+ 1 ) c «+2. (4) 

As can be seen in Fig.|3]also the next-to-leading order coefficient, C4, is strictly positive. 

A new feature shows up in the expansion coefficients at <^(/i 6 ). The coefficient c& is 
positive only below Tq and changes sign in its vicinity. If this pattern persists for higher 
order expansion coefficients one may conclude that the irregular signs of the expansion 
coefficients for T > Tq suggest that the radius of convergence of the Taylor series is not 



related to critical behavior at these temperatures, whereas it determines a critical point 
for T < Tq if also the higher order expansion coefficients stay positive. 

Ratios of subsequent expansion coefficients provide an estimate for the radius of 
convergence of the Taylor expansion, 



The expansion coefficients d n for the quark number susceptibility have been analyzed 
recently for unimproved staggered fermions tl26fl up to n = 6. It has been shown that an 
accurate determination of these expansion coefficients requires large physical volumes. 
Based on a finite volume analysis the radius of convergence has been estimated from 
the Taylor series of the quark number susceptibility to be /ig ~ 180 MeV lEfill . As 
the expansion coefficients d n are directly related to the expansion coefficients c„ of the 
pressure the radius of convergence coincides in the limit n — > °°. For finite n, i.e. n ~ 6, 
estimates based on the Taylor series for the pressure will be about 30% higher than those 
based on the Taylor series for susceptibilities. 

While the quark number susceptibility is expected to diverge at the chiral critical 
point, the pressure, of course, will stay finite. Although % q rises rapidly with increasing 
jJ-q/T, this is partly due to the rapid increase of the pressure itself. A quantity reflecting 
the relative magnitude of fluctuations is given by the ratio, 



This quantity only depends on ratios of Taylor expansion coefficients, and should vanish 
at a 2 nd order phase transition point. It is shown in Fig. |7J As can be seen the ratio rises 
rapidly across Tq and approaches the ideal gas value. It, however, does not show any 
sign of a drop in the vicinity of Tq that could be taken as evidence for the existence of a 
second order transition. For this reason it has been argued in Ref. [ 18] that conclusions 
on the radius of convergence drawn from estimators p n may be premature at least for the 
quark masses used in this calculations. As these masses were significantly larger than 
those used in Ref. 112611 both results are not necessarily in contradiction. Like in the case 
of studies of Lee- Yang zeroes a more detailed analysis of the quark mass dependence 
also is needed in this case. 
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FIGURE 7. The dimensionless ratio &.p/T 2 % q for two values of the chemical potential fi q /T '. Horizon- 
tal lines show the expected results for a hadronic resonance gas and an ideal quark-antiquark gas below 
and above 7b, respectively. 
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